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Abstract 

During cell division, sister chromatids are segregated by the mitotic spindle, a bipolar assembly of 
interdigitating antiparallel polar filaments called microtubules. Establishing a stable overlap region 
is essential for maintenance of bipolarity, but the underlying mechanisms are poorly understood. 
Using a particle-based stochastic model, we hnd that the interplay of motors and passive cross 
linkers can robustly generate partial overlaps between antiparallel hlaments. Our analysis shows 
that motors reduce the overlap in a length-dependent manner, whereas passive cross linkers increase 
it independently of the length. In addition to maintaining structural integrity, passive cross linkers 
can thus also have a dynamic role for size regulation. 

PACS numbers: 87.16.Ka, 87.16.Nn, 87.16.Ln, 87.10.Hk 
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Eukaryotic cells store their chromosomes in a nucleus. During nuclear division, called 
mitosis, the duplicated chromosomes segregate. This process relies on the mitotic spindle, a 
bipolar structure of interdigitating microtubules pp. These hlamentous polymers are phys¬ 
ically cross linked by specihc proteins to maintain the spindle’s integrity. Some of these 
cross linkers are molecular motors that can generate mechanical stresses from the hydrol¬ 
ysis of ATP, which in turn can result in relative sliding between microtubules. Whereas 
the gross architecture of spindles is conserved from yeast to human, their detailed internal 
organisation varies largely between species. 

In spite of its vital importance, the physical principles underlying formation and size 
regulation of mitotic spindles are largely unknown. Some works emphasize the inherent 
ability of microtubules and molecular motors to self-organize into spindle structures by 
means of their mechanical interactions and microtubule assembly [2H1]. Microtubules turn 
over on a time scale of tens of seconds with kinetic differences at their two ends, denoted 
as plus and minus. Other authors have emphasized the role of external factors in spindle 
assembly and maintenance, for example, through macroscopic gradients of various regulatory 
proteins [5HS]. 

On the scale of individual hlaments, notably the regulation of microtubule length by 
length-dependent depolymerization [lOHTT] was studied. Less is known about mechanisms 
regulating the overlap between interdigitating microtubules, which is key for maintaining 
spindle bipolarity in many organisms. Cross linking molecular motors moving towards the 
microtubules’ plus ends tend to reduce the overlap region, see Fig. [^. This might be com¬ 
pensated either by hlament growth at the plus end na or by the action of antagonistic 
motors moving towards the minus end [16]. Both mechanisms, though, turn out to require 
hne tuning of parameters. Passive cross linkers, that is, without motor activity, were long 
thought to merely provide structural integrity and to effectively increase the friction be¬ 
tween sliding microtubules d- In vitro experiments, however, have indicated a regulatory 
function for the passive cross linkers MAP65/PRC1/Asel through recruitment of specihc 
motor proteins [TH| or through generating friction that depends on the overlap size HI. 
Theoretical work showed that these cross linkers and molecular motors can segregate along 
microtubules [20] . Passive cross linkers were also shown to lead to an increase of the overlap 
between hlaments [SI 122]. Consequently, passive cross linkers might play a dynamically 
much more interesting role for spindle homeostasis than previously assumed. 
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FIG. 1: (Color online) Stochastic model for the dynamics of antiparallel microtubules in the 
presence of molecular motors and passive cross linkers, a) Schematic representation of single and 
double bound motors (black circles) and passive cross linkers (red squares), b) Kymograph for the 
motor and cross-linker distributions on one of the filaments with its plus end at z = 1. Empty 
sites are black, sites occupied by motors in green, and sites occupied by passive cross linkers in 
red. Darker colors indicate single bound, lighterr colors double bound particles. The white line 
indicates the position of the other filament’s plus end. c) Distribution of the overlap length for the 
simulation in (b). Simulations are for = 0.316 and 0 ;° = 10“^. For parameter values see text. 

In this work we show that passive cross linkers can oppose the action of molecular motors 
to separate two antiparallel microtubules resulting in the formation of a stable partial over¬ 
lap, see Fig. The overlap length can be tuned by varying the relative concentrations of 
motors and passive cross linkers. A mean field analysis shows that the overlap results from 
a motor-induced sliding stress that decreases with decreasing overlap and an antagonistic 
sliding stress induced by the passive cross linkers. The latter is independent of the overlap 
length and follows from an asymmetry in cross-linker binding at microtubule ends. Finally, 
we investigate the effects of steric interactions between motors and passive cross linkers. 

We describe a pair of antiparallel microtubules by two oriented lattices each with N 
sites, see Fig. Motors and passive cross linkers are each incorporated as particles with 

two identical ‘heads’ that are connected by a linear spring of respective stiffnesses fcm and 
kp. Each site of the lattices can accommodate at most one motor and one passive cross- 
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linker head. Motivated by in vitro experiments, we consider reservoirs of motors and passive 
cross linkers, implying that receptive sites are occupied by particles from the reservoir with 
respective constant rates and 0 ;°. They are occupied with heads of particles already 
bound to the other lattice at respective rates and Wp, where the position of a free head is 
distributed normally around the attached head’s position with width a. The detachment of 
single bound particles occurs at rates and Wp. Assuming slip bonds, the detachment rates 
of a single head of a cross-linked particle increases exponentially with the spring extension 

that is, p exp {fcm^p |^| //m,p}, where /m^p are characteristic force scales for motors and 
passive cross linkers. 

Bound heads can move along a hlament. Motor heads bound to a single hlament hop 
at rate 7 to the neighboring site in the direction of the lattice’s plus end. Heads of passive 
cross linkers hop at rate D to one of the two neighboring sites. Hopping is only possible, if 
the target site is not yet carrying a head of the same species. For cross-linked particles, the 
hopping rates depend on the spring extension ^ as 7 exp and D exp {±fcp,^//p}. 

The sign is chosen such that hops reducing the spring extension are more likely. As for 
the passive cross linker Asel [19], we consider the microtubule ends to constitute diffusive 
barriers: hopping off the lattice is suppressed. For motor particles, the boundary sites are 
treated as bulk sites. 

To determine the motion of the microtubules, we compute the total force generated by 
the extensions of all springs of cross-linking particles. This force is balanced by the friction 
force pn, where v is the hlament’s velocity and /i its friction coefficient. We neglect possible 
Brownian motion of the hlaments. This completes the description of our model. 

We performed stochastic simulations, see Fig. [^,c, for parameters obtained in vitro for 
the motor Eg5 and the passive cross linker Asel: cj^ = 0.0017s“^, ca™ = 1.56s“^, and 
7 = 12s“^ [ini [23l [2l|. For the hopping rate of passive particles we chose D = 12s“^ [29] . 
The cross-linking rates caj^p and the spring stiffnesses kp are not known experimentally. 
Throughout this work, we use Up = 1.2s“^, = 8.4s“\ and fcm,p = O.llpN/nm [30]. For 

the characteristic forces we took the stall force of Eg5, /m,p = 9pN iza, and a is taken to 
be 8 nm. Each lattice had 100 sites corresponding to a microtubule length of 800nm and the 
friction coefficient p was O.lpNs/nm [3T]. The rates w^ p depend on the concentrations of 
proteins in the reservoir and will be used as control parameters. From now on, we scale all 
rates by 7 , lengths by 8 nm, the size of a tubulin dimer, and forces by /m. 
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FIG. 2: (Color online) Overlap length as a function of the occupancy rates and uip from 
stochastic simulations (a) and from the mean-field analysis (b). For stochastic simulations, the 
overlap length has been sampled 10^ times. The crosses indicate the occupancy rates used in 
Fig.[T}3,C, other parameter values are as given in the text. 

For = 0.316 and 0 ;° = 10“^, we observe a stable overlap between the lattices of ap¬ 
proximately 50 sites, see Fig. [^,c. The distribution of overlap lengths is roughly Gaussian 
with a standard deviation of 4 sites. The phase diagram in Fig. shows that for a hxed 
concentration of passive particles, there is an interval of motor concentrations for which a 
stable partial overlap is generated. If their concentration is below a critical value, both hla- 
ments overlap fully. This is in line with the reports in Refs. [211122] that passive cross linkers 
alone maximize the overlap between two hlaments. Remarkably, also for motor concentra¬ 
tions above a critical value, we observe full overlap. This indicates that the ability of motors 
to separate antiparallel microtubules is reduced with increasing motor concentration. These 
observations are reminiscent of the dependence of the particle current in the totally asym¬ 
metric exclusion process (TASEP), which initially increases with the particle concentration 
and then decreases as it approaches a completely hlled lattice |2S]- From Figure]^, one sees 
that the region of partial overlap shrinks with the concentration of passive particles. 

To understand the origin of the stable overlaps, we turned to a mean-held analysis of our 
system. In steady state, the distribution of motors and passive cross linkers along the two 
hlaments are the same when plotted from the plus to the minus end for each hlament. We 
will focus on this symmetric situation and consider the densities Pm,p(3;) of particles bound 
to only one hlament at a distance x from the plus end. Note, that x is now a continuous 
variable. In addition, 0m,p(a^,O denote the densities of cross linking particles, where one 
head is bound at a distance x from the plus end of one hlament and ^ is the spring extension. 
The second head is thus located a distance ^ — x + ^ from the plus end of the other hlament. 
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The time evolution of pni,p is given by the continuity equation 



( 1 ) 


Here, the currents for motors and passive cross linkers are, respectively, given by 


jm — qPm(l Pm Qm) 

jp D iyQpdxpp PpdxQp 5 


( 2 ) 

(3) 


where Qm,p = f Qm,p with = 0jn,p(x, + 0m,p(^ ~ 2 ; + ^, ^) denote the total densities 


of cross linking proteins at position x. The overlap length evolves according to pi = F^ + Fp 
with 


L 



(4) 


0 



exp {—^^/2}/v^^. Here, Pm,p and Qm,p are the corresponding densities on the opposite 
microtubule. 

The densities 0m,p of cross linking proteins also evolve according to the continuity equa¬ 
tion. The expressions for the currents and sources follow the same logic as for Pm,p, but are 
rather involved, such that we refrain from giving them here. 

The mean-field analysis reproduces the essential features of the stochastic simulations, 
see Fig. The densities Pm,p and Qm,p are essentially constant in the overlap and non¬ 
overlap regions of the hlaments [32], see Fig. ii-c. The corresponding force densities are 
constant in the bulk of the overlap region and show pronounced maxima at its boundaries, 
see Fig. |^,d. The total force exerted by cross linking motors is dominated by the bulk 
and tends to decrease the overlap. As the overlap region shrinks, this total motor force 
decreases. In contrast, for passive cross linkers, the force density vanishes in the bulk. The 
corresponding total force is given by the boundary values and independent of the overlap 
length, see Fig. |^, and also tends to increase the overlap. Summing up the forces yields a 
stable steady state overlap length. 

At the level of individual particles, what causes the total net force of the passive cross 
linkers? Consider a passive cross linker in the bulk. The distribution of the spring extensions 
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FIG. 3: Mean field analysis of the steady state. a,c) Densities of single bound (dashed lines) and 
double bound (solid lines) motors (a) and passive cross linkers (c). b,d) Average force density for 
motors (b) and passive cross linkers (d). The dotted lines indicate the end of the overlap region. 
Parameters are as in Fig. [3.,c. 

will be symmetric around zero due to the fast diffusion of bound cross linker heads and 
because the cross linking rate is symmetric with respect to At the boundary, however, 
this distribution is asymmetric because binding sites are absent beyond the hlament’s plus 
end. This leads to a net force. The same effect exists for motors. However, usually the bulk 
forces will mask it. Only in case the motor density is so high that motor movements are 
essentially blocked by steric effects does the end effect dominate. This asymmetry explains 
the maximal overlaps observed for very low and very high motor densities observed in Fig. 

So far, we have not considered steric interactions between motors and passive particles. 
Due to their different sizes, they might indeed occupy different protofilaments on cross 
linked microtubules. Still, a priori, one cannot exclude steric interactions between the two 
protein species. Extending our model to include steric interactions, we hnd that stable 
partial overlaps can be generated under these conditions, too, see Fig. |^. Interestingly, 
for a similar average overlap length, the distribution is broader in the presence of steric 
interactions, see Fig. |^. Up to a critical value of the occupancy rate the phase diagram 
looks similar to the case without interspecies steric interactions, see Fig. Bi. Above that 
value, we observe a phase with fluctuations that are of the order of the filament length, see 
Fig-Hj. More work is necessary to characterize this phase in detail and to identify the origin 
of this behavior. 
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FIG. 4: (Color online) Stochastic model with interspecies steric interactions, a) Kymograph as in 
Fig. & and for the same parameter valnes. b) Kymograph for io^ = 10^ and 0 ;° as in (a). Dashed 
white line indicates the position of the other filament’s minus end. c) Distributions of overlap 
length for (a) (dark grey) and for Fig. (light grey), d) Overlap length as a fnnction of and 
(Up. The grey region indicates highly fluctnating states as (b). The cross indicates the occupancy 
rates used in (a), the dot those used in (b). 

In summary, our analysis shows, that in the presence of motors and passive cross linkers 
stable partial overlaps can be generated between antiparallel hlaments and that the overlap 
size can be tuned by the concentrations of motors and passive cross linkers. Previous work 
has shown that, in the presence of steric interactions between motors, sliding can also be 
induced between parallel hlaments [271 EH]. It will be interesting to explore the consequences 
of these collective effects for the organization of ensembles of microtubules and their possible 
impact on spindle formation. 
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